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1. Introduction 



Almost thirty years after the first serious attempts at formulating a quantum theory of 
gravity [1], still very little is known about what such a theory might be, beyond the semiclassical 
approximation. It is non-renormalizable as a quantum field theory, and there is a general feeling 
that one must seek a non-perturbative formulation of quantum gravity, since the basic concept 
of a smooth geometry becomes incompatible with Hcisenberg's uncertainty principle at very 
small scales. 

Unfortunately, the non-perturbative analysis of nonlinear dynamical systems is a noto- 
riously difficult problem. Quantum gravity is far from the only case where non-perturbative 
effects are important. Other examples include some of the most pressing problems at hand, 
such as climatology and plasma physics. Certainly, the challenge of developping new methods 
for this class of problems must be accepted. Theoretical work in the field has split along two 
currently non-intersecting trajectories; one deals with chaos, the other is related to knots (or 
more generally, things that cannot be made to vanish by a succession of small deformations). 
Although both are studied in the context of gravity, knots appear to be most closely related to 
quantum gravity and the issue of the "small scale structure" of spacetime: relations between 
quantum gravity and knots have appeared both in 3 + 1 dimensions [2], and in the model of 
(2 -|- l)-dimensional gravity [3], which we will pursue in this article. 

Before we begin, it should be stressed that there is still no candidate for a consistent 
theory of "quantum gravity" , nor a proof that such a theory exists, so in a way our problem 
is worse than, say, climatology. On the other hand, numerous authors have suggested that 
(2 -|- l)-dimensional gravity could be used as a simpler model, in which one could formulate a 
quantum theory and investigate issues of interpretation, the choice of time, etc. [4]. 

Witten's work on the subject has led to a complete set of Heisenberg picture observables 
for the quantum theory, an important step in the quantization programme [5]. In spite of this, 
much remains to be done before one could say that quantum gravity in 2+1 dimensions has 
been solved. The authors who have undertaken to follow up Witten's work have encountered 
major obstacles in two different directions. Nelson and Regge, trying to make Witten's reduced 
phase space explicit, found that they could do so up to genus two only by finding the ideal 
of a complicated system of relations among traces, some of which are cousins of the Cayley- 
Hamilton relations for SO (3) matrices [6]. On the other hand, the task of extracting the 
explicit time dependence of the observables (i.e., solving the Heisenberg equations for some 
choice of internal time) was solved in the case of genus one universes by Moncrief [7], and in 
the context of Witten's variables, by Carlip [8]. Interestingly, these authors diffcrred as to what 
the Hamiltonian should be, in particular Carlip's proposed Hamiltonian was unbounded from 
below. Following Moncrief 's remarks to this effect, Carlip showed that his formulation led to 
a "Dirac square root" of the Whceler-DeWitt equation, when one demands invariance of wave 
functions with respect to the mapping class group [9] . 

Our purpose in this article is to pursue the quantization of (2 + l)-dimensional gravity to 
the point of providing a complete set of Schroedinger observables, a Hamiltonian, a consistent 
choice of ordering of the operators and, in computable form, wave functions for various types 
of scattering problems. We hope that this will help to formulate some of the fundamental 
problems of quantum gravity in terms of specific, computable questions. We stress that some 
of these "problems" involve the very consistency of the formalism which we are proposing, so 
we do not claim to have proved that the theory exists. 



Many authors have worked on (2 + l)-diniensional quantum gravity, and it sometimes 
seems that each one has a different approach. We will briefly list some of these approaches in 
terms of two main schools of thought. The issue is how to reduce the infinite dimensional phase 
space by exploiting the infinite number of gauge symmetries, to obtain a finite-dimensional 
system which could be quantized exactly. One approach is to fix the gauge by choosing a 
slicing of spacetime in constant curvature slices (York gauge [10]), or some other choice of 
slicing - we will call this the ADM school. The other relies on the fact that (2 + l)-dimensional 
spacetimes can be regarded as a particular sheet of solutions (with "maximal Euler class" ) of 
the Chern- Simons gauge theory for an ISO(2, 1) connection - the ISO(2, 1) school. 

1- The ADM School 

1.1 Particles and cones. 

The geometry of an K^ universe with particles is a multi-cone, i.e. a flat surface with 
conical singularities, where the deficit angle at each singularity is related to the mass of the 

particle. The two-particle classical dynamics was solved explicitly by Deser, Jackiw and 
't Hooft [11], the quantum scattering of a light particle by a heavy one was worked out by 
Deser and Jackiw [12], and the general two-particle quantum scattering problem was solved by 
't Hooft [13]. The three-particle classical problem was solved by Lancaster and Sasakura [14]. 

1.2 Riemann surfaces. 

The geometry is a genus g Riemann surface, from a constant extrinsic curvature foliation 
of spacetime. The surface coordinates are chosen so that the metric is conformally flat (Hosoya 
and Nakao) or conformally hyperbolic (Moncrief). To compute the Hamiltonian, one solves an 
equation for the conformal factor. In the first case, it is not a coordinate scalar so the equation 
must be solved on coordinate patches - it is not known whether it admits a solution for genus 
g > I. In the conformally hyperbolic case, Moncrief finds a Lichnerowitz equation for the scalar 
conformal factor, and proves the existence of a solution for all g [7]. The Shroedinger equation 
is easily written down for genus one [15], and the solutions can in principle be found from the 
modular invariant eigenfunctions of the Laplacian on TcichmuUcr space, which are known as 
the Maass forms [16]. 

The virtue of the ADM approach is that it is most similar to (3 -I- l)-dimensional gravity. 
Its drawback is that it encounters some technical problems that are most similar to those of 
(3 + l)-dimensional gravity. 

2- The IS0(2,1) School 

2.1 Reduced phase space quantization. 

One finds a complete set of phase space observables which commute with the constraints. 
The reduced phase space is related to the moduli space of flat ISO{2, 1) connections (Witten, 
[5]). The problems begin when one attempts to parametrize this space and compute the Poisson 
brackets: one can form an infinite set of ISO{2, 1) invariants from the holonomy, but these 
are not independent variables: They are related by non-linear constraints, some of which are 
cousins of the Cayley Hamilton identities. Nelson, Regge, Urrutia and Zertuche [17] have 
succeded in identifying the reduced phase space explicitly up to genus 5 = 2. The canonical 
quantization process leads, in some cases, to quantum groups [18]. The reduced phase space 
can be interpreted as a Hamilton-Jacobi formulation [19], so the representations of the reduced 
algebra should give the wave functions in the Hcisenberg picture. No attempt is made to 
compute the Hamiltonian or to formulate quantum dynamics questions. 



2.2 Covariant quantization. 

One works with the IS0(2, 1) "homotopies" along loops which form an arbitrarily chosen 
basis of TTi{T,g), and later imposes invariance with respect to a change of basis, or mapping 
class transformation (Carlip, [20]). The difficulties encountered by Regge et al. can be avoided 
by using these homotopies rather than the invariants derived from them, leaving out a global 
ISO{2, 1) symmetry. The homotopies can be parametrized in the "polygon representation" of 
(2 + l)-dimensional gravity, and the symplectic structure is known (H.W., [21]). Mapping class 
invariant scattering amplitudes can be computed by the method of images, by summing over all 
mapping class images of the 'in' state. Two problems are encountered: Each term in the sum 
can only be computed if one knows the Hamiltonian, for a given choice of time. The maximal 
slicing choice only allows one to compute the Hamiltonian explicitly only for genus g — I (see 
ADM school). The other problem is that the sum is over the non-abelian, infinite, mapping 
class group, and it is not clear a priori which terms might be small, or if the sum converges. 



2.3 Planar Multiple Polygons 

Recently, 't Hooft has proposed another approach based on a slightly different polygon 
representation, where one chooses a piecewise flat surface, which leads to a number of planar 
polygons, with 5*0(2, 1) identifications of the boundary edges. There is a remarkably simple 
symplectic structure and Hamiltonian for this system; since these variables are directly related 
to a choice of generators of the first homotopy group, the mapping class group acts non-trivially 
on the wave functions [22]. 

This article can now be summarized in one line: We combine our earlier calculation of the 
Hamiltonian in the polygon representation, with the method of images and the "stationnary 
phase theorem" , to advance along the covariant quantization programme [23] . 

The reader is referred to the previous article on the polygon representation of (2 + 1)- 
dimensional gravity for the canonical variables, the Hamiltonian, and most importantly, the 
relation between the polygon representation and the ISO{2, 1) homotopies [19]. We will begin 
this article directly with the canonical quantization in Sec. 2. The mapping class group in- 
variance is discussed in Sec. 3, where we write down invariant amplitudes for various types of 
scattering problems. In Sec. 4 we give the conclusions and discuss the problem of time. 



2. Quantization in the Polygon Representation 

In this section, we will draw heavily on results from the classical theory; their derivations 
can be found in [19] and references therein, and will not be repeated here. The reduced phase 
space of (2 + l)-dimensional gravity can be parametrized hy 2g + N three-vectors E{^) and 
as many 5*0(2,1) matrices M{fi), with the following Poisson brackets and constraints (/i — 
l,2,..,2g + N). 

{E-{^,),E\^,)}^e'^'^E,{^i) (2.1) 

{E-{^l),M\{^i)}^e'^'>''Mac{p) (2.2) 

{M\{^Ji),M-^{^,)}^Q (2.3) 

The vectors E°-{ii) together with their identified partners —M~^{ii)E{ii) form a closed polygon, 

J = Y,iI-M-\fi))E{fi)^0 (2.4) 

and the cycle condition for the SO {2, 1) identification matrices leads to the constraints 

P- = ie-^-Wefc « (2.5) 

where 



W = (m{1)M-^{2)M-^{1)M{2U ■ ■ ■ ( M(2g - l)M-^{2g) 
M-^{2g - l)M{2g)\M{2g + 1)M(25 + 2) • • • M{2g + N) 



(2.6) 



The constraints J°' and P" generate a Poincare algebra with the brackets (2.1)-(2.3). The 
masses of the N particles are given by the constraints 

iJ(/i) = P^i^i) + sin^ (n{fi)) w (2.7) 

which generate translations of each particle along its worldline. 

It is convenient to define the canonical variables P"(Af) = (l/2)e"''^Mc6(^) and 

x,.,.^(.,.,.,,.,-aMd^^M) 

where 



Ji^i) = {I - M-\fi))E{fi) (2.9) 

The variables X{^) and P{fJ.) have canonical Poisson brackets: 

{Paifl), X'ifl)} - Sa' (2.10) 

The constraints can be written in terms of the canonical variables, using J{ij-) = X{fi) A 
P(/i) and 

M\ = S\ + P,e-\ + (Jl + P^ - l\ (5\ - ^^^ (2.11) 

for /i = 1, ..., 2(7 (for which M{p) is hyperbolic (a boost) [24]), or for /i = 2g + 1, ..., 2g + N, 

M\ - S\ + P,e-\ + (Vl - P^ - l) (s\ - ^^\ (2.12) 

The SO {2, 1) constraints become 

J=^X(/.)AP(A*)«0 (2.13) 

while the translation constraints P w are defined implicitly in terms of P{fi) by (2.11) and 
(2.12). The explicit form of the function P(P{fi)) w can be derived from the relations (2.11)- 
(2.12). 

We will assume the universe is a closed surface with genus g > I; for the other cases 
the arguments below carry over substitutes for the "time" and "Hamiltonian" , the appropriate 
expressions drawn from the previous article. We propose the following choice of "internal time" , 
in a frame where M(l) is a boost in the {y i)-plane: 

T= K (2.14 

P^(l) ^ ' 

The Hamiltonian is given in terms of the variablcsM(/i), /i > 1, by solving the constraint 

W = ( M(l)M-i(2)M-i(l)M(2) j • • • ( M(2.g - l)M-^{2g)M-^{2g - 1)M(25) j « / (2.15) 

for tr{M{l)), where M(l) is a boost of magnitude b with axis x. One finds 

P,(2)Q,-P,(2)Q, 
^= P,^(2)-P,^(2) (2-1^) 

where 



1 



= — e 



abc 



M(2)(M(3)M"i(4)M"i(3)M(4)) ■ 



2 

(2.17) 

(M(2g - l)M-\2g)M-\2g - l)M{2g)) ) 

/ cb 

There are two residual "translation" constraints and one "boost" constraint. The former 
come from M{1)M-^{2)M-^{1) « (M(2)M(3)M-i(4)M-i(3)M(4) • • • M{2g + N))-^, 

P"^(2)-Q^«0 (2.18) 

P^(2)-Q^w0 (2.19) 

The boost generator is 



J-^^ /-M-i(Ai) E'^{fj,)^0 (2.20) 

The remaining phase space variables are E{p,), M{p,) for /i = 2, ..., 2g + N. Their time evolution 
is generated by the Hamiltonian (2.16) and any linear combination of the first-class constraints 
(2.18)-(2.20), which commute with the Hamiltonian. 

If TV 7^ 0, one can choose a gauge where the vector X{2g + 1) is aligned with the a:;-axis 
and X{2g + 2) vanishes at all times; M{2g + 2) is then a pure rotation with angle ^{2g + 2), 
and 

T^^:p±^ (2.21) 

P^(25 + l) ^ ' 



P"" =€'"'"'{ M{2g + 3)---M{2g + N)M{l)---M{2g)M{2g+l)] wO (2.22) 

cb 



py EEey'"'fM{2g + 3)---M{2g + N)M{l)---M{2g)M{2g+l)\ «0 (2.23) 



H = Tr (M(2g + 2) • • • M{2g + N)M{l)M-^{2) ■ ■ ■ M-^{2g - l)M(2.g) j (2.24) 

Note that the internal times (2.14) and (2.21) are linearly related, as well as any choice 
which is linear in the E's, since cPE{^i)/dt — 0. Such times are also "inertial", in the sense 
that the position of a particle as seen from the observer, E{^), follows a straight trajectory at 
constant "velocity". 
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We begin by constructing the quantum theory in the most naive way, then pick up the 
problems as they appear (see the mapping class group, next section). We quantize the canonical 
phase space X{^), P{ij) in the usual way {^ — 2, 3, ..., 2g + N). 

P^i^) - -^ (2.25) 

[Pa(M), X\^l)] = -tSa' (2.26) 

To prove the equivalence of this quantization to one which could be carried out in the 
variables E{ii), M{fi), one must express these variables in terms of X(/i) and P{p), and show 
that there exists an ordering of these expressions such that the commutators form a represen- 
tation of the algebra (2.1)-(2.3). An acceptable ordering is obtained by placing all the X's on 
the right: 

E{^,) = P MP AX)- ^^-^P{P ■ X) (2.27) 

where the index /i ~ 2, ...,2g + N has been omitted in the right hand side expressions, and 
M^/^, the "square root" of M{ii), was given explicitly in the previous article, Eqns. (6.7) and 
(6.9). 

As mentioned, the ordering of the operators in (2.27) leads to the correct commutator al- 
gebra for E{fx),M{fj,). For example, using P"(^) = ~ir]''''d/dX''{p.), one finds [P"(Ai), E''{fi)] = 
—ie°''"^Ec(p). This solution to the ordering problem is not unique. 

The Wheeler DeWitt quantization would consist in taking the set of square integrable 
functions on IR ^^ ^ ', which are annihilated by the constraint operators: 

^^^i{Xi^l)■,^l^l,■■■,2g + N}) (2.28) 

J"V = (2.29) 



P"V = (2.30) 

Equations (2.29) state that V' must be a Lorentz scalar in the variables X{^). The other 
constraints are the Wheeler- DeWitt equations. These can in principle be obtained explicitly 
from (2.11), (2.12) by replacing everywhere Paip) by the operator —ihd/dX°'{ii). The result 
is a non-rational function of the derivatives, not likely to give a physically acceptable quantum 
theory. For this reason, we will prefer to partly fix the gauge as indicated above, by choosing the 
internal time (2.21) and the corresponding Hamiltonian (2.24), and quantize in the Schroedinger 
picture: 

'^^ - H{{P{^^)■,^^ - 2, • • • , 2g + N})^ (2.31) 

One must still impose the two "momentum constraints" (2.22)-(2.23), 



^ -1 

^xbc I 



e''""'iM{2g + 3)---M{2g + N)M{l)---M{2g)M{2g+l)\ iP ^ (2.32) 

ey'"'(M{2g + 3)---M{2g + N)M{l)---M{2g)M{2g+iyj V = (2.33) 

and require that the wave function be a Lorentz scalar. 

The solution to the Schroedinger equation, for an initial state tlj{X{^), 0), /x = 2, 3, ..., 2^ + 
N , is given by 

^{{X{^,)}■ T) = e^THi{-^a/^xi^)} ^({X(/,)}; 0) (2.34) 

The Haniiltonian operator is defined by the Fourier decomposition: 

^TH{{-^^/^X{^L)} ^iKo-X ^ ^TH(Ko) ^2.35) 

Finding initial states which satisfy the constraints is not difficult in principle. Only the 
constraints (2.32) and (2.33), which depend on the momenta in a non-polynomial way, appear 
to present a difficulty, but in P-space these become algebraic relations, which summarize the 
cycle conditions, for which explicit solutions can be constructed (this is not always easy). 



3. Mapping Class Group and Quantum Scattering 



The preliminary theory of (2 + l)-diniensional quantum gravity which we constructed in 
Sec. 2, is formally similar to the non-relativistic quantum mechanics of a set of free particles 
propagating in a flat three-dimensional background, but with an unusual "kinetic energy" 
function. This non-interacting picture can be understood as follows: We have defined the 
polygon variables by cutting up and unfolding the manifold, to get a polygon embedded in 
Minkowski space. In this classical picture, the corners of the polygon follow straight timelike 
lines - so indeed they behave very much like free particles! 



3.1 The Mapping Class Group 

The variables which we are using are associated to a particular choice of generators of the 
fundamental group for a genus g surface with N punctures. A different choice of generators can 
lead to a different set of variables for the same spacetime. The group of transformations from 
one set of generators to another, called "mapping class group", reflects a discrete symmetry 
of the polygon representation. One must demand that the wave function be mapping-class 
invariant, which leads to interesting "scattering" effects [23]. 

Before we begin to construct invariant wave functions, we will review the mapping class 
group and its action on the polygon variables, including "time" and the Hamiltonian. We must 
also make sure that our choice of operator ordering is consistent with this symmetry group. 

The mapping class group is presented by the following generators [25] (the loops ai^hi 
are a standard basis, where each ai intersects only the corresponding bi at only one point; 
Tui , Ty. , Tzi are Dehn twists [20] ) . 



Tui ■ bi — > biUi (3.1) 



Ty : Ui^ Uib^ ^ (3.2) 



{ai -^ aib~^ai+ibi+ia:[j^^ (3.3) 

bi -^ ai+ib~^^a~^-i^biai+ibi+ial^-^ (3.4) 

For genus g = 2, the representation of these generators in ISO{2, 1) was given in [19], and 
can be recovered with the change in notation Ui —^ ai,Vi —^ bi. For example, 

'M(l) (/-M(l))£;(l)-M(l)£;(2) 
P(a0-| Hl-Mil))OA I (3^g^ 

1 
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/M-i(2) (/-M-i(2)-M-i(l))^(l) \ 

p(5,) = +(^ - M-H2))Ei2) + (/ - M-\2))OA ^^^^^ 

so that we find, applying Equation (3.1), that p{ai) is unchanged and p{bi) becomes 



Pibi) 



/M-i(2)M(l) {I ~ M-^{1) - M-\2)M{1))E{1) \ 
+ {I - M-^{2) - M-^(2)M{l))E(2) 
+ {I - M-^{2)M{1))0A 

V 1 / 



(3.76) 



By comparing this to the expression (3.7a), one can deduce the action of the mapping class 
transformation on the variables E(p) and M{p,) (this is not easy). 

r E{1) ^ E{1) + E{2) (3.8) 

I £;(2)^M-i(l)£;(2) (3.9) 

^"' • ] M(l) -> M(l) (3.10) 

Im(2)^M-1(1)M(2) (3.11) 

and similarly for /i = 3 and /i = 4, and for any other element of the mapping class group. 
These relations can also be derived from the polygon representation of the variables E{fj,) and 
M{p) (Figure 3.1). The loop ai is a loop starting at O, going through the segment E{1) at the 
point X, through M^^(1)_E(1) at X' and back to O. Likewise 6i starts at O, goes through the 
segment M~'^(2)_E(2) at F', then continues through Y and returns to O. Now one can deform 
the combined loop 6iOi smoothly into the loop which is represented in Figure 3.1, where a 
small semicircle surrounds the corner B of the polygon. One can then choose a new set of 
cuts which avoids the corner B: The first cut is now E{1) + E{2) (dashed line), while the 
second cut is as before but in a frame which has been transported around the loop E{2) once: 
£■(2) -^ M^^{\)E{2). Any mapping class group transformation can be obtained in a similar 
fashion from the polygon representation, by sliding one corner of the polygon around a loop (a 
Dehn twist). 

Mapping class transformations leave the cycle conditions invariant, since they are auto- 
morphisms of 7ri(S). This implies that the Hamiltonian, which was computed in the previous 
section by solving the cycle conditions in a frame where M{1) is a pure x-boost, is invariant as 
long as the frame is adjusted so that this condition is preserved. We will thus choose the frame 
independently for each term in the sum over the mapping class group. This requires adjusting 
the initial and final conditions appropriately, by an overall 5*0(2, 1) transformation. Time, 
defined as the variable canonically conjugate to this Hamiltonian, is mapping-class invariant up 
to a constant. 

One must check that the mapping class group action is consistent with the choice of 
operator ordering. If one chooses to place the E^s and X's on the right, the bracket algebra, 
as well as the relations that give X{ji) as a function of E{fj,), are preserved under the mapping 
class group transformations (this can be checked explicitly, and derives from the fact that all the 
relations are linear and homogeneous in E and X). This differs significantly from the situation 
in York's extrinsic gauge, where the task of finding an operator ordering consistent with the 
mapping class group generally requires that one consider multi-valued wave functions, with a 
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modular weight which depends on the ordering chosen [9]. We feel that this requirement is 
due to the choice of slicing, and is not a fundamental property of (2 + l)-dimensional gravity. 
Indeed, the polygon variables do not exclude the possibility of choosing the extrinsic curvature 
slicing; rather they do not specify any foliation: We have decoupled the slicing (gauge) from 
the global variables (observable), and simply dropped the former. The difficulty discovered by 
Carlip, is that it is difficult to find a consistent operator ordering for the relations which give 
the ADM operators in terms of the quantized IS0{2, 1) homotopies. 



3.2 Invariant wave functions - The method of images 

One must find mapping class invariant wave functions which solve the Schroedinger equa- 
tion. We first review what has already been achieved in both the ADM and the ISO{2, 1) 
schools. In the ADM approach [7, 26], the Hilbert space is the space of modular invariant 
square integrable functions over TeichmuUer space. To find the solutions to the Schroedinger 
equation, one considers the eigenfunctions of the Hamiltonian operator. The Hamiltonian is 
not yet known for genus greater than one (one must solve a Lichnerowitz equation), but it is 
expected that if a solution could be found, then the ordering problems would be intractable. 
This means that one is probably limited to the case of the torus, where the Hamiltonian is 
known and efforts are under way to compute the wave functions explicitly. 

In the ISO{2, 1) approach, one computes the amplitude for the scattering N particles on 
a plane as a sum of amplitudes over all possible braidings of the worldlines of the particles; 
if one is dealing with a compact universe with punctures, then one sums over the mapping 
class group, rather than the braid group. Each term in this sum is a path integral with the 
Chern-Simons action, over all Poincarc connections on three-manifolds with fixed "initial" and 
"final" conditions, and a given topology and mapping (or "braiding", for particles). Besides 
the difficulty of computing the path integral (which is related to the difficulty of finding the 
Hamiltonian), one must deal with an infinite sum over a non-abclian group; it is not known 
whether the sum can be ordered as a converging series, or if the result is computable. 

In both approaches, one can separate the task at hand into two sub-tasks: 

(1) The dynamical problem (finding the Hamiltonian, or computing a sum over histories with 
fixed mapping). 

(2) The mapping class problem (finding modular- invariant eigenstates of the Hamiltonian 
operator, or performing the sum over the mapping class group). 

Both the dynamical and the mapping class problems have been solved for the scattering 
of two particles on IR [13], and the Hamiltonian is known for the quantum torus [7], with work 
under way on the mapping class problem [33]. 

So far in this article, we have solved the dynamical problem in the general case of a genus 
g surface with N punctures. In what follows, we will attempt to provide a solution to the 
mapping class problem, by constructing the invariant wave functions as a series of computable 
terms. We will give a procedure to decide on which terms (or elements of the mapping class 
group) give a significant contribution to the amplitude, and argue that the number of such 
terms is finite. 

We will consider the amplitude to go from an initial state | Xi >, which we choose to be 
an eigenstate of the "position" operators at Ti — 0, {X'^(p); fi — 2,- ■ ■ ,2g + N}, to a final state 
I X2 > at T2 > 0. Neither of these states is mapping class invariant, but one can construct an 
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invariant states | X >inv by summing over all mapping class images of X: 

\X>^nv- J2 \pi9^)^> (3-12) 

ie{m.c.g} 

This eigncstatc has infinite norm, equal to the cardinality of the mapping class group. We 
will use the formal expression (3.12) as a starting point, and derive physically acceptable wave 
functions below. 



3.3 Scattering Amplitudes 

Even though the spacetimes considered here are not asymptotically flat, one can define 
non-interacting "in" and "out" states as follows. Since there is no potential energy term in 
the Schroedinger equation, the system is non-interacting as long as the mapping class group 
symmetry is not taken into account. It is easy to see from the construction of invariant states 
in (3.12), how the mapping class group affects this non-interacting picture: A realistic state 
might be a superposition of position eigenstates (3.12), with some variance a^. Thus, the 
support of each term in the sum is a fuzzy patch centered at p{gi)X. If the patches for various 
mapping class images overlap, then one expects to pick up interference terms - this is the 
essence of "interactions" in (2 -I- l)-dimensional gravity. We will define "in" and "out" states 
by demanding that such overlaps do not occur - practically, this requires that the X{p) be 
"sufficiently large", so that for every element gi of the mapping class group, the image of X{p) 
under gi lies in the tail of the wave packet, which is exponentially suppressed: 

'i{^i,^}^\X{^Ji)-p{g,)X{^Ji)f»a^ (3.13) 

This criterion depends on the variance of the wave packet; one can get a more universal condition 
if one assumes that the variance is of order fi . In units h — 1 this gives 

V{a., 0, II X{p) - p{g,)X{pL) f» 1 (3.14) 

We are now ready to calculate the scattering amplitude: Taking two asymptotic invariant 
states (3.12) and the Hamiltonian (2.24), and recalling that the mapping class group leaves H 
invariant and translates T by a constant, the amplitude is just 

< 1|2 >- Y. < P{9,){X2,T^)W"^^^-^'^^p{9^){XuT,) > (3 15) 

i,je{m.c.g.} 

^ J2 <X2,r2|e-«(^-^r')|p(g^.-i)p(5,)(Xi,ri)> (g^g) 

i,je{m.c.g.} 

where {X2, T^} — p{gj){X2, T'2}, etc.. Finally, since no loop of 7ri(E) is privileged with respect 
to any other, the products {g~ gi) will cover the mapping class group uniformily, and absorbing 
the infinite factor equal to the cardinality of the group, 

<1|2>^ ^ <X2,T2|e-*^(^^-^i')|p(g.)(Xi,ri)> (3.17) 

ie{7n.c.g.} 
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Thus, the ampUtude at {X2JT2} is a sum over the contributions of each mapping class 
image of the source, {XI, Tl}, propagated to {X2,T2} by the "free-particle" Hamiltonian 
H{P{ij), ^ — 2, ..,2g + N). Note that the sum is infinite, and the expression (3.17) is not 
necessarily well-defined. We will argue below that (3.17) is computable as a distribution on a 
certain limited set of wave functions tp{{Xi,Ti}); it would be very interesting to identify an 
appropriate space of test functions for which the distribution (3.17) is well-defined, but this lies 
beyond the scope of this article - the difficulties arise because the mapping class group does not 
act properly discontinuously on the space of polygons. 

3.4 The stationary phase theorem, and the sum over mapping class images 

We will consider an initial wave packet centered at a suitable Xi, where by "suitable" we 
mean that the images {XI, Tl} have no accumulation point (3e : Vi, X]u(^(m) ~ p{gi)^{p))'^ > 
e). The orientation of the "beam" is defined in a frame-independent way, to get an 50(2, 1)- 
invariant wave function. Wc further require that the wave packet a(K) have support only over 
P's which correspond to solutions of the cycle conditions which form a faithful representation 
of the fundamental group [5], [21] (other sheets of solutions include totally collapsed handles, or 
curvature singularities with a surplus angle equal to a multiple of At:). 

V'(Xi,Ti)= f dKa{K)e'^<^^-^°'^ (3.18) 

where 

i^ = {P"(A*);M = 2,3,---,2g-fiV} (3.19) 

dK == n^ d^P{fi) (3.20) 

The wave function at (X2,T2) is given by propagating the initial state (3.18) with the 
invariant propagator (3.17). With the notation ujk — H{K) and the definition 

{X(,),r(,)} = p{g,){Xi,n) - {X2,T2) (3.21) 

the wave function becomes 

V'(X2,T2)- Y^ /"dii'a(i^)e^(-^-^w-'^'=^w) (3.22) 

ie{ni.c.g.'\ 

The stationnary phase argument implies that the only contributions to the integral are 
from values of K where the phase varies slowly as compared to a{K)] this picks out the trajec- 
tories near the classical solutions 

where the first factor in the r.h.s. is the group velocity of the packet, which can be calculated 
explicitly given the Hamiltonian (2.24). In the scattering problem, the initial state, final state 
and i are fixed, so that (3.23) can be solved for K\ let us denote this solution by K{i). The 
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only elements of the mapping class group which give a large contribution to the sum (3.22) 
are such that K{i) is not in the exponentially suppressed tail of the wave packet: a{K{i)) > e 
for some appropriately chosen e. Since we have assumed that the images of the source have 
no accumulation point, one expects that only a finite number of terms will contribute and the 
wave function (3.22) has a finite norm. This is necessarily a rough statement at this point; a 
general proof of computability of the expression (3.22) for some appropriate Hilbert space of 
wave functions is not yet available. We have examined a number of specific examples to various 
degrees of completeness, and found that the expression (3.22) is computable for initial data 
peaked about sufficiently non-singular polygons. One encounters the following situations. 

(1) If the images are widely spaced and their emissions do not interfere, only one of the classical 
trajectories K{i) lies within the wave packet a{K); one recovers the free propagator. This 
occurs when the images X{i) are sufficiently distant from each other as compared to the distance 
between X2 and Xi, so that the images K{j) different from K{i) do not lie within the support 
of the wave packet a{K). There is no scattering in this case. 

(2) Two classical trajectories connect the initial and final states, such as with geodesies on a 
cone (two-particle scattering). If the initial wave packet is wide enough, e.g. includes both 
P(l) and M(2)P(1), then the scattering amplitude involves the interference between these two 
terms, as in a two slit experiment. Other terms in the sum, such as (M(2))"P(1) for n > 1, 
correspond to the events where one particle winds around the other n times - the corresponding 
K may lie within the packet a{K) for small enough n, but these trajectories are not "nearly 
classical" and will contribute little, by the stationnary phase agument. According to 't Hooft's 
more complete analysis [13], such terms are suppressed by a factor of order (e/27r)'*", where e 
is the deficit angle at the conical singularity. 

(3) One of the braid generators of the mapping class group (winding around a small loop) 
acts as a small translation, so that one picks up interference from a large number of images 
spaced on a line: This amplitude is related to scattering of light by a grating, and describes 
the quantum dynamics of a thin wormhole. If one of the mapping class generators ga modifies 
K{i) very little, where K{i) is a classical solution which lies within the packet a(k), then by 
applying this generator again and again one obtains a sequence of contributions from K(^j^-^, 
with (7(j„) — {ga) "5(*)- The series (3.22) may be exactly summable in this case; in the example 
of the torus the result should be related to the Maass forms [16]. 

(4) Various mapping class generators contribute significantly to the sum: In the general case 
where more than one generator has "almost non-proper" action on K, one must consider the 
interference of the amplitudes propagating freely from a large number of images, which are 
distributed in a complicated fashion. This corresponds to a multiple scattering amplitude, and 
the expression (3.22) is related to the scattering by a complicated multidimensional "grating". 



3.5 Quantum gravity near the "big bang", and through the "big bounce"? 

Besides the scattering from asymptotically free states, another problem of interest is the 
(2+ l)-dimensional equivalent of the "big bang" . Many solutions of the classical theory expand 
from an initial singularity (although smoothly, unlike in 3 -I- 1 dimensions). For A^ = and 
genus g > 2, all solutions have either an initial or a final singularity [27]. This brings up a very 
interesting situation: As one approaches the singularity, the topological features are increasingly 
tiny and a wave packet of given variance will contain an increasing amount of mapping class 
images. This means that the number of terms in the invariant wave function (3.22) increases, 
and it becomes less and less reasonable to interpret the result in the semiclassical picture, by 
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referring to a background geometry with interference terms acting as effective interactions. The 
"geometry" becomes completely fuzzy, but the wave function can still be calculated - this gives 
one a window on the non-pcrturbative regime, and small scale structure, of a theory of quantum 
geometry. 

The wave function (3.22) is probably computable up to any desired level of accuracy at a 
fixed point away from the singularity, although the number of terms which must be considered 
in the sum can be large if one is close to the initial singularity. Short of actually performing this 
calculation, one can draw a few tentative conclusions. The increasing density of mapping class 
images as one approaches the big bang indicates that the inevitable anisotropy of classical uni- 
verses with non-trivial topology is likely to be smoothed out: Any effect which occurs before the 
quantum gravity effects lose significance is likely to appear isotropic to a future observer, even 
in a topologically non-trivial cosmology. By the same argument, it appears plausible that the 
wave function at {X, T} away from the big bang would depend very little on the precise initial 
condition on ^(Xi,Ti). This would be very convenient, since it eliminates the problem of hav- 
ing to specify an initial condition for the very early universe; it places all the burden of selecting 
among the possible histories, on the measurements which are presumably performed later on, 
and collapse the wave function to something resembling a classical background geometry. 

The last compact universe problem which we wish to discuss is the "big bounce" . The 
initial state is a wave packet centered about a contracting universe, which would classically 
contract to a singularity and then expand again on the other side. The semiclassical picture 
breaks down as one approaches the singularity, since the wave function becomes a superposition 
of an increasingly large number of mapping class images and the wave function extends to cover 
all possible universes with the given topology. It is impossible to think in terms of a "geometry" 
in this region, yet interference patterns from the mapping class images is likely to create "dark 
regions" and "bright regions" in the space of geometries. It appears that the large number of 
mapping class images may lead to randomization of the wave packet, and a very non-localized 
wave function on the other side, as in the previous case. This may be related to the chaotic 
behavior of particles close to the singularity [28] . 



3.6 Scattering of Particles and Wormholes on IR 

If the universe has the topology IR with g handles and N punctures, the amplitudes are 
calculated in much the same way but the Hamiltonian function is simpler. The geometry at 
infinityis that of a cone with a helical shift, and the axis of the cone defines the direction of 
the "time" vector. The deficit angle of the conical geometry gives the total energy, H, and the 
helical shift corresponds to the total angular momentum [29] [11]. The cycle condition for closed 
topologies is replaced by a single equation, which gives the overall holonomy for a loop which 
goes around the circle at infinity; we will consider only wave functions with support limited by 
the condition that this overall holonomy is a rotation, since otherwise closed timelike curves 
occur in the semiclassical region (this implies that there is at least one particle, since wormholes 
are always tachyonic [21]). The angle of this rotation is the Hamiltonian: 



H = Arcsin{P) (3.24) 

where 
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pa ^ e"'"' ( M{l)M~^{2)M-^(l)M{2) ■ ■ ■ M{2g + NU (3.25) 

2 V J cb 

This function has been written down expUcitly as a function of the momenta for two 
particles [11] and three particles [14], and can be calculated, with an increasing amount of 
work, for any number of particles. Note that there is an ambiguity in the definition of the 
Hamiltonian, which is only given modulo 2tt. This ambiguity is related to the fact that the 
holonomies are representations of the fundamental group in the covering group of S0{2, 1). It 
is usual to restrict one's attention to the case where the holonomy M{ij) associated to each 
individual particle (/i = 2g + l,...,2g + N) lies on the component simply connected to the 
identity, and similarly for the overall holonomy at infinity - this corresponds to particles of 
reasonably small mass [21], and so-called "physically reasonable" conditions at infinity [30]; the 
consequence of these asumptions is that for the g — Q case, one can show that there exists 
a global spacelike slicing (no time machines) [31]; this fact was first stated without proof in 
[11]. We will enforce these conditions by choosing wave packets with support limited to these 
"reasonable geometries", for g = Q. In the general case, there is no avoiding this ambiguity, 
and one must decide if the wave function is multivalued or if it transforms non-trivially under 
a 27r rotation [23]. As before, we will make the choice of a scalar wave function and postulate 

*(X,T)- ^ /dii:a(i4:)e*(^-^(*)^""=^c)) (3.26) 

i^{m.c.g.^ 

This expression is similar to Carlip's proposal for exact scattering amplitudes [23], the 
only difference being that we are considering wave packets, rather than the pure propagator. 
This allows us to appeal to the stationnary phase theorem and argue that the sum is limited 
to a finite number of mapping class images, which are such that the classical trajectory which 
connects an image to {X,T} corresponds to a momentum K(i) which lies within the support 
of the packet a{K). It would be very interesting to see an explicit solution of the three-particle 
scattering problem in a specific situation, by using (3.28) and the expanded expression for the 
Hamiltonian [14]. 
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4. Conclusion; the "problem of time' 



In this article, we have developped the covariant quantization programme which was 
initiated by Witten and Carhp, using the "polygon representation" of the reduced phase space. 
In this representation, the classical (2 + l)-dimensional gravity appears formally similar to a 
set of free particles propagating in H^^^. We gave the explicit expression for the Hamiltonian 
function and quantized the theory canonically. The mapping-class invariant wave function 
can be written as a sum of freely propagating amplitudes, where each term represents the 
propagation from one mapping-class image of the initial state. Interactions occur when these 
amplitudes interfere. This can produce interference patterns analogous to either the two-slit 
experiment, or the diffraction by a regular grating, in some special cases. 

It is not obvious that the sum over mapping class images is computable for wave functions 
which have support over singular geometries. Since singular geometries occur rather generically 
in (2 + l)-dimensional gravity, it would be of great interest to prove that the sum is well-defined 
for a specific Hilbert space of wave functions. 

Also of interest, is the opportunity to compare this quantum theory for the torus {g = 1), 
to other author's results in the extrinsic time. Joining our results and those of Carlip [9], we have 
indirectly shown that the theories are closely related if one considers only static issues, since 
it was shown that the ADM, Witten-Carlip and Polygon variables are different representations 
of the same reduced phase space, and that the relations between these representations are 
consistent with the operator orderings chosen [9]. However, it is not clear whether the quantum 
dynamical theories, given by Schroedinger's equation in York's time or in our internal time 
T, are equivalent. We will argue below that they are not, but first we must define what we 
mean by "equivalent" . We will use the term "observable" in the traditional sense of quantum 
mechanics, without requiring that they should be constants of the motion. An example of an 
observable might be, e.g., L^ = E'^(p). 

Definition Two quantum systems, described by the wave functions ip and ip' in times t and 
t' , are said to be "equivalent", if for any set of observables Ai, {i — 1, ...,N), the expectation 
values < Ai >t and < Ai >t' follow the same trajectories in IR . 

Let {Ci, H} be a complete set of commuting observables, where H is the Hamiltonian for time 
t, and let K be the Hamiltonian for time t'. The wave functions ip and ip' can be expanded 
in a basis of eigenstates of the commuting observables. We now explain why the theories of 
(2 + l)-dimensional gravity in different times are not equivalent. 

(1) K must be diagonal in this basis (if it were not, then the expectation values of the observables 
in the state ip' would not be constant and the equivalence criterion would be violated). 

(2) The initial state ip'{t' = 0) must have the same probability for each CSCO eigenstate as 
tp{t — 0) (a trivial consequence of the definition above: set d = Ai). 

(3) The wave functions tp(t) and ip'{t') differ only by the phase factor of each coefficient in 
the expansion in the basis of eigenstates. The time-dependence of these phases in the two 
formulations are given by (pnif) = 'P'ni^) + ^nt' and 0„(t) = 0„(O) + Hnt. One easily checks, 
by comparing the form of the Hamiltonian for our choice of internal time T to that which 
corresponds to the extrinsic time, that these phases do not follow the same trajectory in IR , 
where d is the number of commuting observables. This further implies that one can easily find 
a set of A^ = 2 observables (which do not commute with H) such that their expectation values 
do not follow the same trajectory in IR . Thus, either the quantum theories in different times 
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give different dynamics, or they cannot be compared in the sense that it is impossible to set 
up the same experiment in both pictures. The attempts at avoiding this fundamental problem 
by appealing to the path integral formulation of quantum mechanics, as a unified formalism 
to which both Schroedinger quantizations are presumably equivalent, usually end up hiding 
the problem in either an inadequate definition of the boundary conditions, or of the measure: 
as Guven and Ryan have emphasized [34], two inequivalent theories cannot be simultaneously 
equivalent to a third one! 

Until one can gain a better grasp on this problem, and find an appropriate Hilbert space 
for which the wave function (3.22) is computable, it is impossible to claim that this canonical 
formalism for (2 + l)-dimensional gravity is acceptable as such; probably it will be necessary 
to look first at the simplest case g = 1, where it is possible to compare the results with those 
of the ADM school. 

One might consider extending this work to 3 + 1 dimensions along the lines of the exactly 
solvable theories proposed by Horowitz [35]. A reduced phase space for Horowitz's topological 
B/\F theory can be found [36, 37], and it is probable that the canonical quantization programme 
can be carried out along the same lines as those laid out in this article. 

The canonical quantization programme of (2 + l)-dimensional gravity was faced with 
two problems: that of computing the Hamiltonian, without which quantum dynamical issues 
cannot be addressed, and that of finding wave functions that are invariant under the action of 
the mapping class group. The first problem has essentially been solved, although this required 
abandoning the extrinsic gauge which is favored by the ADM school for its relevance in 3 + 1 
dimensions; as we argued above, the change of internal time is not a trivial step in a quantum 
theory! As for the second problem, an explicit construction of the invariant wave functions 
was given. Even if the computability cannot be established in the general case, we hope that 
this article has provided a practical scheme with which one can tackle specific non-perturbative 
problems of quantum gravity in 2 + 1 dimensions. 
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Figure Captions 

Figure 3.1 The loop aibi crosses twice through the edges of the original polygon, at X = X' 
and at y — y . By recutting the surface in a different way prior to unwrapping it as a polygon, 
one obtains a first cut which goes directly from A to C (dotted line), so that the loop aibi 
crosses only once an edge of the new polygon, exactly in the same way as oi crossed one edge 
of the original polygon. The switch from one set of cuts to another, and from ai to aibi, is an 
example of a mapping class transformation. 
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